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Abstract. In this paper, we consider the compressible Euler-Maxwell equations arising in semi- 
conductor physics, which take the form of Euler equations for the conservation laws of mass density 
and current density for electrons, coupled to Maxwell's equations for self-consistent electromagnetic 
field. We study the global well-posedness in critical spaces and the limit to zero of some physical 
parameters in the scaled Euler-Maxwell equations. More precisely, using high- and low-frequency 
decomposition methods, we first construct uniform (global) classical solutions (around constant equi- 
librium) to the Cauchy problem of Euler-Maxwell equations in Chemin-Lerner's spaces with critical 
regularity. Furthermore, based on Aubin-Lions compactness lemma, it is justified that the (scaled) 
classical solutions converge globally in time to the solutions of compressible Euler-Poisson equations 
in the process of non-relativistic limit and to that of drift-diffusion equations under the relaxation 
limit or the combined non-relativistic and relaxation limits. 

Key words. Euler-Maxwell equations, classical solutions, Chemin-Lerner's spaces, non-relativistic 
limit, relaxation limit 

AMS subject classifications. 35L45, 76N15, 35B25 

1. Introduction and main results. The increasing demand on semiconductor 
devices has led to the necessity of a deep and detailed understanding on the mathe- 
matical theory of various charge-carrier transport models. Of these important mod- 
els, the classical hydrodynamic model (also named as the Euler-Poisson equations), 
which treats the propagation of electrons in semiconductor devices as the flow of a 
compressible charged fluid in an electric field, has received increasing attention. For 
the cases of high electric field and submicronic devices, the Euler-Poisson equations 
of fluid dynamical form can represent a reasonable comprise between physical accu- 
racy and reduction of computational cost in real applications, the reader is referred 
to [22] for more explanations. When semiconductor devices are operated under some 
high frequency conditions (such as photoconductive switches, electro-optics, semicon- 
ductor lasers and high-speed computers), magnetic fields are generated by moving 
electrons inside devices, then the electrons transport interacts with the propagating 
electromagnetic waves. In this case, the transport process is typically governed by the 
Euler-Maxwell equations, which is more accurate than the Euler-Poisson equations, 
since the electromagnetic field obeys Maxwell's equations instead of Poisson equation 
for the electric field only. 

After some appropriate re-scaling, the compressible Euler-Maxwell equations are 
written, in nondimensional form, as 



for {t, x) e [0, +oo) X R^(A^ > 2). Here the unknowns n,u— {ui,U2, ■ ■ •, un)^ , E = 
{Ei,E2, ■ ■ •, i?Ar)^, B = (_Bi, i?2, • • •, -Bjv)^(T transpose) denote the electron density. 
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dtn + W ■ {nu) ^ 0, 

dt{nu) + V • (nu (g) u) + VP(n) 



n(E + £u X B) - sii, 



(1.1) 



< eA^9tE — V X B = enu, 
edtB + V X E = 0, 
A^V-E^n-n, V • B 0, 
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electron velocity, electric field and magnetic field, respectively. The pressure P{n) 
satisfies the usual 7-law 

(1.2) =Fo7i^(7> 1), 

where Pq > is some physical constant. The system (II. ip is called isentropic if 
7 > 1 and isothermal if 7 = 1. r, A > are the (scaled) constants for the momentum- 
relaxation time and the Debye length, c — (eoi^o)"^ > is the speed of light, where eo 
and vo are the vacuum permittivity and permeability. Setting e — The independent 
parameters t, A and e which arise from nondimensionalization, are assumed to be very 
small compared to the reference physical size. The symbols V, •, x and (g) are the 
gradient operator, the scalar products, the vector products and the tensor products 
of two vectors, respectively, n > is the doping profile, which stands for the density 
of positively charged background ions. 

It is not difficult to see that the above Euler-Maxwell equations consist of a 
quasi-linear hyperbolic system, the main feature of which is the finite time blowup 
of classical solutions even when the initial data are smooth and small. Hence, the 
qualitative study and device simulation of (jl.ip are far to be trivial. In this paper, 
our main aim is to establish the global well-posedness and justify some singular limits 
for the Cauchy problem. For this purpose, the Euler-Maxwell equations are 
equipped with the following initial conditions for n, u, E and B: 

(1.3) (n, u, E, B)(x, 0) = (no, Uq, Eq, Bo)(a;), x € R^, 
which satisfies the compatible conditions 

(1.4) A^V • Eo = n - riQ, V • Bq = 0, xe R^. 

1.1. Singular limit analysis. It is convenient to state previous works and main 
results of this paper, we first introduce some singular limits in the scaled Euler- 
Maxwell equations at the formal level, including the non-relativistic limit, relaxation 
limit as well as combined non-relativistic and relaxation limits. 

Firstly, we observe the non-relativistic limit {i.e. e ^ 0). Let r = 1 = A and 
(n^, u*^, E"^, B"^) be the solution of the following equations 

dtn^ + V • (n^u^) = 0, 

= ~n%W + eu" X B^) - n^u", 
^"^■^^ I edtW - V X B^ = er^'^u^ 

edtB^ +V xW = 0, 
V-W^n-n^, V-B^ = 0. 

Formally, we see that the limits u", E° of n^, u"^, E^ as e — > satisfy 

(1.6) \ at(nOuO)-|-V- (nOuO(guO)-l-VP(n") = -nOEO-nOuO, 

[ V • E° = fi - n", V X E° = 0, 

which is the well-known Euler-Poisson equations for semiconductors. The irrotational- 
ity of E° implies the existence of a potential function such that E*' — —V <1>'~'. Then 
using the Green's formulation, (jl.6p can be reduced to the form of the conservation 
law with a non-local source term, e.g., see |13j . 
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Secondly, we justify the relaxation limit (i.e. r 0) in the Euler-Maxwell equa- 
tions The diffusion limit was first introduced by Marcati and Natalini [3T] for 
the Euler-Poisson equations (II. 6p . Set e = 1 = A. To do this, as in [21], we define the 
following scaled transform 

(1.7) in\ u^ B-)(<, x) = (n, ^u, E, b) . 

Then the new variable (n"^, u"^, E^, B^) satisfies 
dtn^ + V ■ (n^u^) = 0, 

T^dtin^'u^) + r^V • (n^u^ (g) u^) + VP{n^) 

= -n^{W + TU^ X B^) - n'vJ, 

TdtB-" -f V X E^ = 0, 



Formally, the limits A/", £ of n'^, E^ as r 
tions 



satisfy the so-called drift-diffusion equa- 



(1.9) 



dtN^V ■ {VP{N)+N£), 
V ■£ = n~M, Vx£ = 0, 
N{Q,x) = no, 



which is a system of diffusion equations for the electron density, and maintains the 
parabolic-elliptic character. 

Lastly, we study the combined non-relativistic and relaxation limits in the Euler- 
Maxwell equations (|l.ip (i.e. £,r — > 0). Set A = 1. From the "©(l/r) time scale" in 
(|1.7p . where the superscript r is replaced by (r, e), the new variable (n''^'^-', u'^'^''^\ e'^'^\ 
B^^'^)) satisfies 



9tn(^'^) + V • (n(^'^)u(^'=)) = 0, 

r2at(n(^'=)u(^'^)) + r^V • {n'^'^'^^u'^^'^^ (g) u^^'^)) + VP(n(^''^)) 



(1.10) <^ 



TedtB'^^'^^ + V X E^^'"^) = 0, 

V • E(^'^) ^n- , V • B(^'=) = 0. 



j(r,e)^(r,e) 



Obviously, in the process of combined limits r, e 
also satisfy the drift-diffusion equations (|1.9p . 



0, the limits TV, £ of n(^'=), E 



1.2. Main results. In the past ten years, the Euler-Poisson equations (|1.6p have 
attracted much attention. There are many contributions in mathematical analysis, 
such as the well-posedness of steady-state solutions, global existence of classical or 
entropy weak solutions, large time behavior of classical solutions, relaxation limit 
problems and so on, the reader is referred to [11 11 [iTJ [12l [13l [Ml [15l [19l [21] and 
the references therein, also including ourselves [lOl [27l [28] , while the Euler-Maxwell 
equations are much more intricate than the Euler-Poisson equations, not only because 
of Maxwell's equations, but also because of the complicated coupling of the Lorentz 
force (E + u X B). In contrast, not so many works have been devoted to the study of 
Euler-Maxwell equations. Up to now, only partial results are available. 
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Using the Godunov scheme with the fractional step and the compensated com- 
pactness theory, Chen, Jerome and Wang [5] constructed the existence of a global 
weak solution to the initial boundary value problem for arbitrarily large initial data 
in L°°(R). In [16], assuming initial data in Sobolev spaces i/*(R'^) with higher regu- 
larity (s > 5/2), a local existence theory of smooth solutions for the Cauchy problem 
of non-isentropic Euler-Maxwell equations, where the pressure-density function ()1.2p 
is replaced with the energy equation, was established by modificating the classical 
semigroup-resolvent approach of Kato [T7] . In [231 [Ml 121] , based on the existence 
theory of Kato and Majda [I71[5Dj, Peng and Wang justified the non-relativistic limit 
(e — > 0), the quasi- neutral limit (A — >■ 0) and the combined non-relativistic and 
quasi-neutral limits (e = A 0) for the Euler-Maxwell equations in virtue of 

the analysis of asymptotic expansions. Their results show that the Euler-Maxwell 
equations converge towards the Euler-Poisson equations, e-MHD system and incom- 
pressible Euler equations in some time-interval independent of the parameters e and 
A, respectively. 

However, the well-posedness and singular limits for the Euler-Maxwell equations 
in several dimensions are still far from well-known, in particular, in the frame- 
work of critical spaces. In the present paper, we shall answer this problem. More 
concretely speaking, we shall consider a small perturbation near the constant equilib- 
rium state (n, 0, 0, B) which is a particular solution of the Cauchy problem (|l.ll) - (|1.3p . 
and obtain the global existence and uniqueness of classical solutions. We choose the 
critical Besov spaces in space-variable x as the basic functional setting, where the 
regularity index (cr = 1 + is just the limit case of classical existence theory of 

Kato and Majda [17l[2^. Although this idea has been used to study the compressible 
Euler-Poisson equations (|1.6I) in [TUl [571 [2E] recently, it should be pointed out that 
the Euler-Maxwell equations are essentially different from (|1.6p . In comparison with 
the methods in [TOl [27l [28] , we have to face with several technical difficulties arising 
in the uniform a priori estimates of classical solutions in critical spaces. The first one 
is lack of the low- frequency estimate of magnetic field B, which does not lead to the 
exponential decay near equilibrium in view of the standard definition of norm of Besov 
spaces. Another one is that the nonlinear terms (pressure, Lorentz field, etc.) will 
hinder us establishing the uniform estimates with respect to the singular parameter 
couple (t, e). To overcome these difficulties, we add the new content in the proof of 
the local existence and (uniform) global existence of classical solutions. Actually, the 
Chemin-Lcrner's spaces L^(i3p^) in 4^ are introduced, which is a refinement of the 
usual spaces L^{B^ ^), and some uniform frequency- localization estimates in Chemin- 
Lerner's spaces with critical regularity are established, for details, see Lemmas l3.3ll3.4l 
and Lemmas l3.6ll3.7l Based on the uniform estimates, we further rigorously justify the 
singular limit problems for ()l.l|) - (jl.3|) in Sect. Il.ll bv the standard weak convergence 
methods and the application of compactness theorem in [26^. 

Throughout this paper, the regularity index a ^ 1 + N/2. First of all, we state a 
local existence and uniqueness theorem of classical solutions to the Cauchy problem 
(|l.ip - p.3l) away from the vacuum. 

Theorem 1.1. Let n > be a constant reference density and B G be any 
given constant. Suppose that UQ — n, Uq, Eq and Bq — B £ i(R^) satisfy no > and 
the compatible conditions (T^, then there exist a time Tq > and a unique solution 
(n, u, E, B) of the system j|) - l77S)) such that 



(n,u,E,B) e C^([0,To] X R^) with n > for all t e [0,To] 
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and 

{n - n, u, E, B - B) e CT^Bh) n C^^lS^^'). 



Remark 1.1. To avoid excessive commutators arising from the nonlinear pres- 
sure term by using Fourier frequency-localization method, we introduce a function 
transform in Sect. 13.11 such that the Euler-Maxwell equations (|l.ip is reduced to a 
symmetric hyperbolic system. Based on the previous effort in |10| . we obtain the 
local existence of classical solutions in the Chemin-Lerner's space with critical regu- 
larity (Proposition 13.11) . Theorem 11.11 follows from Proposition 13.11 and Remark 13.11 
readily. As a matter of fact, the new result is applicable to generally symmetrizable 
hyperbolic systems, which enriches and develops the classical existence theory of Kato 
and Majda [ITlE^. 

In small amplitude regime, we get the uniform global well-posedness of classical 
solutions to the Cauchy problem ()l.ip - (|1.3p in critical spaces. From now on, we set 
the scaled Debye length to be one (A = 1). 

Theorem 1.2. Let n > be a constant reference density and B e be 
any given constant. Suppose that no — n, uo,Eo and Bq — B G ^^(R^) satisfy 
the compatible conditions There exists a positive constant Sq independent of 

singular parameter couple (t, e), such that if 

\\{no - n,Uo,Eo,Bo - B)||bj^ < (Sq, 

then there exists a unique global solution (n, u,E,B) of the system ifj. j|) - l77]g)) satis- 
fying 



(n,u,E,B) G C^([0,oo) x R 



N\ 



and 

{n - n,u,E,B - B) e C{Bl,{'R.'')) nC\B^^j\K'')). 
Moreover, the uniform energy estimate holds: 
n,u,E,B - B) 



(1.11) < ColKno - n, uo, Eq, Bq - B)||bj^ 









) 


+ 











for < T, £ < 1, where the positive constants fJ-o,Co are independent of (t, e). 

Remark 1.2. Together with Theorem ll.il Theorem 11.21 directlv follows from the 
standard continuation argument and the crucial energy estimate (jl.lip which presents 
the dissipation rates of all the components in the solution. Noticing that the coupled 
electromagnetic field (E, B) appears in the nonlinear source terms of Euler system, 
which indeed does not affect the character of corresponding linearized form, so we 
can take the full advantage of "Shizuta-Kawashima" skew-symmetry condition which 
was well developed for general hyperbolic systems of balance laws [18j [29] to capture 
the dissipation rate of density function, see Lemma 13.41 On the other hand, from 
the proof of Lemmas 13.6113.71 we see that the electromagnetic field generated by the 
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compressible electron flow exhibits a weak dissipation property, which is essentially 
different from the pure Maxwell's equations, although the low-frequency estimate of 
dissipation rate of B is absent. In addition, we track the singular parameters t and e 
in the proof of (jl.lip . which plays a key role in the study of related limit problems. 

As a direct consequence of Theorem II. 2 1 we can obtain the large-time asymptotic 
behavior of global solutions near the equilibrium [fi, 0, 0, B) in some Besov spaces. 

Corollary 1.3. Let (n, u,E,B) he the solution in Theorem ll.B. it holds that 

\\n{-,t) - n,u(-,t),E(-,t)||^„_e'(j^„^ ^ 0, \\B{-,t) -B\\g^-r-,' ^^^^ 0, 

as the time variable t — > +oo, where p = {N > 2) and e' > 0. 

Remark 1.3. Recalling the proof of Corollary 5.1 in [9], Corollary [T3] is followed 
by a minor revision. The definite convergence rate to the equilibrium (n, 0, 0, B) will 
be studied in the future work. 

Next, we state the non-relativistic limit of uniform global solutions to (|l.ip - (|1.3p 
for any fixed momentum relaxation time r > 0. 

Theorem 1.4 (Non-relativistic limit). Let t — 1 and (n*^, u*^, E"^, B^) be the 
global solution of iLl\) - [T7^) given by Theorem lLSi Then there exists some function 
(n°,u*',E'') which is a global solution to the Euler-Poisson equations il.6\) satisfying 
(n° - n,uO,E°) G C([0, oo), i(R^)) such that as e->-0, it holds that 

(r^^u^ VeE'^) ^ (n°,u",0) strongly m C{[0,T], {B^r^\R^))ioc), 

VB'^^O strongly in L'^{B^^(R^)), 

(E^ B^ ) ^ (E" , B) weakly* in L^{Bl^{R^)), 
for any T > and S G (0, 1). Moreover, it yields 

||(r.°-n,u",E°)(t,.)bj^(R«) 

(1.12) < Ci||(no-rl,uo,Eo,Bo-B)||Bj^(R«), i > 0, 

where Ci > is a uniform constant independent of e. 

Secondly, we justify the relaxation limit for the Euler-Maxwell equations (|l.ip . 
To this end, we consider the Cauchy problem for the re-scaled system (|1.8|) subject 
to the initial data 

(1.13) {n\u\ W, W){Q, x) = (no, ^uq, Eq, Bq) {x). 

It follows from Theorem 11.21 and the "©(l/r) time scale" (|1.7p that there exists a 
unique global in-time classical solution (n"^, u^, E'^, B'^) to the system (|1.8[) and (|1.13[) . 
Then, we have 

Theorem 1.5 (Relaxation limit). Let e = 1 and {n^ ,W ,W) be the global 
solution of i f_?.($|) and ULIS]) obtained from Theorem \1.2l Then, there exists a function 
{Af,U,£) which is a global solution to the drift- diffusion equations ll.9\) satisfying 

[NM, £) e C([0, oo), B^i(R^)) X l2([o, oo), B£i(R^)) x C([0, oo), B^i(R^)) 
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such that as T ^ 0, it holds that 

(n^,T2u^, V7E^) (7V,0,0) strongly m C([0, T], (B^_7'^(R^))ioc), 

weakly in L^{Bl^(R^)), 

VB^^O strongly m L|(B^_7^(R^)), 

(E^,B^) ^ (5,B) weakly* in L??(B^_i(R^)), 
for any T > and 6 e (0, 1). Moreover, it yields 

||(AA-n,£)(i,.)||B|.,(R") 
(1-14) <C2||K-fi,uo,Eo,Bo-B)||sj^(K«), t>0, 

where C2 > is a uniform constant independent of t . 

Finally, what left is the combined non-relativistic and relaxation limits for (jl.ip . 
From Theorem [L2] and ()1.13p where the superscript t is replaced by (t, e), it is shown 
that there exists a unique global in-time classical solution (rS'^''^\ u^'^'^^ E'^'^''^\ B'^'^''^^) 
to the system (|1.10l) and (I1.13p . Furthermore, we get 

Theorem 1.6 (Combined non-relativistic and relaxation limits). Let 
(,i(T,£)^u(T,e)^ j3(r,e)^B(T,e)) fee the global solution of fZjD)] and ITJ^) obtained from 
Theorem \l/A Then, there exists a function {J\f,U,£) which is a global solution to the 
drift- diffusion equations U.9\) satisfying 

(AA,W,£) GC([0,^),i?2%(R^)) X l2([0,oo),B2-,i(R^)) xC([0,oo),i?f,i(R^)) 

such that as T — > and e ~> simultaneously, it holds that 

(n("'^),r2u("'^),V7?E("^"))^ (AA,0,0) strongly tn C{[0,T], {B^\r''))ioc), 

u'^^'^^^U weakly m i|(B^^i(R^)), 

VB(^'^)^0 strongly in L^B^^\li^)), 

(E(^.^),B(^''^)) (£,5) weakly* in L??(B^_i(R^)), 

for any T > and S G (0, 1). Moreover, it yields 

\\{J\f-n,£){t,-)\\B^jRN) 
(1.15) <C3||(no-n,uo,Eo,Bo-B)||sj^(R«), t>0, 

where C3 > is a uniform constant independent of (r, e) . 

Remark 1.4. To the best of our knowledge, these limit results (Theorems II. 4llL6|) 
show the convergence globally in time, which have not been appeared in the published 
literatures. In comparison with that in f^^llMlllS], they hold true in the functional 
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spaces with relatively lower regularity, which can be regarded as a supplement to the 
theory of singular limits for the Euler-Maxwell equations Let us mention that 

the combined non-relativistic and relaxation limits obtained in Theorem II .61 does not 
require any (communication) restriction between r and s. That is, one can fix any of 
the two parameters t and e and let the other tends to zero, which is the genuinely 
combined limits. 

Remark 1.5. It is worth noting that Chemin-Lerner's spaces are first introduced to 
establish the uniform a priori estimates with respect to (r, e) and justify the combined 
limits. As a matter fact, this approach developed by the current paper can be applied 
to study other limit problems with two (or more) independent singular parameters. 

Remark 1.6. There is no additional conceptual difficulty in considering the tem- 
perature effects and the corresponding balance equation {i.e. non-isentropic Euler- 
Maxwell equations), although the estimates are quite tedious. 

The rest of this paper unfolds as follows. In Sect. [51 we introduce the Littlewood- 
Paley decomposition and recall the definitions and some useful results on Besov spaces 
and Chemin-Lerner's spaces. Sect. |3]is devoted to the proofs of main results, which 
is divided into five subsections for clarity. In Sect. 13.11 we first rewrite the Euler- 
Maxwell equations (|1.1|) as a symmetric hyperbolic system in order to obtain the 
effective a priori estimate by using Fourier frequency localization. Furthermore, we 
give the local existence of classical solutions in Chemin-Lerner's spaces with critical 
regularity. Then in Sect. 13.21 we deduce a new uniform a priori estimate under 
some smallness assumption, which is used to achieve the (uniform) global existence 
of classical solutions. Sect. 13. 3[ Sect. 13.41 and Sect. 13.51 are in turn dedicated to the 
justification of the non-relativistic limit, relaxation limit as well as combined non- 
relativistic and relaxation limits of Euler-Maxwell equations. 

Notations. Throughout the paper, C stands for a uniform positive constant 
with respect to (t, e). The notation f ~ g means that / < Cg and g < Cf. Denote 
by C{[0,T],X) (resp., C^{[0,T],X)) the space of continuous (resp., continuously dif- 
ferentiable) functions on [0,T] with values in a Banach space X. For simplicity, the 
notation \\{f,g,h,k)\\x means + + ||^||x + ||fc||x, where f,g,h,k e X. We 

omit the space dependence, since all functional spaces (in x) are considered in R^. 
Moreover, the integral fdx is labeled as J / without any ambiguity. 

2. Tools. The proofs of most of the results presented in this paper require a 
dyadic decomposition of Fourier variable. Let us recall briefly the Littlewood-Paley 
decomposition theory and the characterization of Besov spaces and Chemin-Lerner's 
spaces, see for instance [Sj [7] for details. 

Let ((/?, x) be a couple of smooth functions valued in [0, 1] such that (p is supported 
in the shell C(0, |, |) = e R^|| < |^| < §}, x is supported in the bah S(0, |) = 
UeR^liei < |}and 

q=0 

Let S' be the dual space of the Schwartz class S. For f £ S' , the nonhomogeneous 
dyadic blocks are defined as follows 

A^if:=x{D)f = h*f with h = T-'x, 



AJ ■.^^{2-m)f = 2''^ hi2'^y)fix~y)dy with = V, if 9 > 0. 
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Here *, J-^^ represent the convolution operator and the inverse Fourier transform, 
respectively. Note that h,h £ S. The nonhomogeneous Littlewood-Paley decomposi- 
tion is 

f=J2\f in S'. 

q>-l 

Define the low frequency cut-off by 

p<q-l 

According to the above Littlewood-Paley decomposition, thus we introduce the ex- 
plicit definition of Besov spaces. 

Definition 2.1. Let 1 < p < oo and s e R. For 1 < r < oo, the Besov spaces 
Bp J, are defined by 

5](2^1|A,/|U.r)'<cx) 

«>-l 

and Bp ^ are defined by 

I^Bl^^ sup 2«^||AJ||iP <^. 

q>-l 

Let us point out that the definition of B^ ^ does not depend on the choice of the 
Littlewood-Paley decomposition. Now, we state some classical conclusions, which 
will be used in subsequent analysis. The first one is Bernstein's inequality. 

Lemma 2.2. Let fc e N and < Ri < i?2. There exists a constant C, depending 
only on i?2 and N , such that for all 1 < a < b < oo and J £ L"^ , we have 

Supp J-/ci5(0,i?iA)^ sup <C'=+iA'=+^(i-*)||/|U.; 

\a.\ — k 

Supp J-/cC(0,i?iA,i?2A)^C-'=-iA'^||/||L« < sup ||6)"/|Ua <C^+ia'=||/||l», 

\a\ = k 

where T j represents the Fourier transform on f . 

The second one is a compactness result for Besov spaces. 

Proposition 2.3. Let 1 < p,r < oo, seR and e > 0. For all (f) € , the 
map f ^ (j)f is compact from B'^'^/ to B^ ,^. . 

On the other hand, the study of non-stationary partial differential equations re- 
quires spaces of type Lf^{X) := Lp{0,T;X) for appropriate Banach spaces X. In 
our case, X is expected to be a Besov space, so the fundamental idea is to localize 
the equations through the Littlewood-Paley decomposition. Then it is easy to obtain 
L^(LP) estimates for each dyadic block. Performing a (weighted) £^ summation is the 
most natural next step. But, in doing so, we get bounds in spaces which are not type 
Lf^{Bp^) (except if p = r). This leads to the definition of Chemin-Lerner's spaces 
first introduced by J.-Y. Chemin and N. Lerner ^4 , which is the refinement of the 
spaces Lip{Bp ^). 

Definition 2.4. For T > 0, s G R, 1 < r, p < oo, set (with the usual convention 
if r — oo) 

q>-l 
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Then we define the space L'^{B^ ^) as the completion of S over (0,r) x by the 
above norm. 

Furthermore, we define 

Ct{b;j :=Z??(B;,jnc([o,T],B;j 

and 

cUb;,,) := {/ e c\[o,nB;^,)\dj e Z??(s;,,)}. 

The index T will be omitted when T — +00. Let us emphasize that 
Remark 2.1. According to Minkowski's inequality, it holds that 

II./IIZp (B.^) < '^ir>p, II/IIzp (B.^) > '^ir<p. 

JThen, we state the property of continuity for product in Chemin-Lerner's spaces 
Proposition 2.5. The following estimate holds: 

II-^5|IZp(B= J < ^(II/IIl?,1(L-)II.9IIi-2(b._^) + \\9\\Li;?iL^)\\f\\lPf^B; J 

whenever s > 0,1 < p < cxd, 1 < P, Pi, P2, Ps, Pi l£ 00 and 

11 11 1 

P Pi P2 P3 Pi 

As a direct corollary, one has 
whenever s > N In, - = — + —. 

— I p pi ' P2 ^ 

In addition, the estimates of commutators in Li^{Bp ^) spaces are also frequently 
used in the subsequent analysis. The indices s,p behave just as in the stationary case 
[3 llOj whereas the time exponent p behaves according to Holder inequality. 

Lemma 2.6. Let 1 < p < 00 and 1 < p < 00, then the following inequalities are 
true: 

2««||[/,A,]Ag||L^(L.) 

[ C'c,||/||zpi(b.+1)II.9|Izp2(b.^), s = n/p, 

where the commutator [■, ■] is defined by [f,g] ~ fg — gf, the operator A = div or V, C 
is a harmless constant, and Cq denotes a sequence such that ||(cq)|j/i < 1, j; = + 

Finally, we state a continuity result for compositions (see [T|) to end up this 
section. 

Proposition 2.7. Let s > 0, 1 < p,r, p < 00, F £ wl^l^^''^{L;R) with F(0) = 
0,T e (0, 00] and v G L^{B^^r) ^ L^{L°°)- Then 
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3. The proofs of main results. In what follows, we focus on the proofs of 
main results. For clarity, we divide them into several subsections, since the proofs are 
a bit longer. 

3.1. Reformulation and local existence. In this section, we reformulate 
(|l.ll) - (|1.3p in order to obtain the effective a priori estimates by means of Fourier 
frequency localization. 

For the isentropic case (7 > 1), let us introduce the sound speed 



(3.1) g=^^[^{n)-^), 



B 



and set -0 = ili{n) corresponding to the sound speed at a background density n. Define 

2 

Set 

W:= (e,u,E,F)T. 

Then the system (|l.ip can be reduced to the symmetric hyperbolic system for smooth 
solutions: 

dtQ + V'divu = — u • V Q — ^^^^gdivu, 

dtn + i^VQ+'^ = -u- Vu- ^pVp- (E + eu x (F + B)), 



(3.2) 



9tE - iV X F = nu + h{Q)n, 
StF + |v X E = 0, 



V • E = -/;,(£i), V-F = 0, 
1 _i — 2 

where h{Q) = {(Po7)^^ (^^t* + 'i/;)}T-i — is a smooth function on the domain 
{q\-^^Q + "0 > 0} satisfying /i(0) = 0. The initial data (11.31) becomes into 

(3.3) I^lt^o = (^'o,uo,Eo,Fo) 

with 

2 

7- 



eo-^-Y(V'("o)-V^), Fo = Bo-B. 



Under the symmetrization transform (13. ip . the initial data (13.31) satisfies the corre- 
sponding compatible conditions 



N 



(3.4) V • Eo = -/i(£'o), V-Fo = 0, x e R 

Remark 3.1. The variable change is from the open set {(n,u,E,B) e (0,+oo) x 
R^ X R^ X R^} to the open set {W^ e R x R^ x R^ x R^|3^g + > 0}. It is 
easy to show that for classical solutions (n, u, E, B) away from vacuum, (ll.ll) - (|1.3p is 
equivalent to dS^-dSS]) with ^^g + -4}>Q. 

For the isothermal case where 7 = 1, the form of (|3.2I) is still valid with ip = \/Pq, 
while the symmetrization transform depends on the following enthalpy variable change 

(3.5) g = ■\/i^(lnn — Inn), 



for the details, see e.g. [TU]. 
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Without loss of generality, we shall study the system p.2p - p.3p for 7 > 1 and 
prove main results, since the case of 7 = 1 can be discussed in the same way. 

In [T71 [50], Kato and Majda established a local existence theory for generally 
symmetric hyperbolic systems pertaining to data in the Sobolev spaces with higher 
regularity. Recently, using the regularized means and compactness argument, we have 
established a local existence in the framework critical Besov spaces for the Euler- 
Poisson equations ()1.6|) . see |10) . In the present paper, we further strengthen the 
result such that it holds in Chemin-Lerner's spaces with critical regularity. Our result 
reads as follows. 



7-1 
2 



Proposition 3.1. For any fixed < r, e < 1, assume that Wq £ 1 satisfying 
go + Tp > and ^3.4\ l. then there exist a time Tq > ( depending only on the initial 
data) and a unique solution W to ll3.2\) - \3.3]) such that W G C^([0,ro] x R^) with 
^g + ip> for all t G [0,ro] and W e Ct„{BI^) n C^^{B^X^). 

Proof. Let the assumptions of Proposition 13.11 be fulfilled. By a proper revision, 
the local result in [T^ can also be adapted to the Euler-Maxwell equations (|3.2l) - p.3p . 
That is, there exist a time Tq > (depending only on the initial data) and a unique 
solution W to (I32I)-(IS31) such that W G C^i[0,To] x R^) with ^g + i/i > for all 
t G [0,To] and W G C([0, Tq], 1) n C^P, Tq], Bj,! In order to prove Proposition 
01 it suffices to show that W G £??(Bf.i) and Wt G Z^f (B^J^). 

Indeed, applying the operator Aq{q > —1) to (|3.2p . we infer that {Aqg, A^u, A^E, 
AqF) satisfies 



(3.6) 



dtAqg + V'A^divu = -(u • V)Ag(? + [u. A,] • - ^Aq{gdWu), 
dtAqU + i)AqVg + iAgU = -(u • V)AgU + [u. A,] 

-^Aq{gVg) - AgE - eA^u x B - £Aq{u x F) 
dtAqE - iV X A,F = fiAqXi + Aq{h{g)vi), 
dtAqF + iV X A,E - 0, 



where the commutator [•, •] is defined by [/, 5] ~ fg ^ gf ■ 

Then multiplying the first equation of Eq. p.6p by Aqg, the second one by A^u, 
and adding the resulting equations together, after integrating it over R^, we have 
the energy equality 

i|(||A,e||i.. + ||A,u||i.) + i||A,u||i. 

= i y" divu(|A,^)|2 + |A,un + j ([u, A,] • ^gAqg + [u. A,] • VuA,u) 

+ j Aqg{\/g- Aqu)-^^ j {[Aq, g]AwnAqg + [Aq, g\W g ■ AqXi 



(3.7) - J AqE-AqU-e J Aq{uxF)-AqU, 



where we have used the fact e(A^u x B) • AqU — 0. 

On the other hand, multiplying the third equation of Eq. p.6p by -^AqE and 
the last one by 4AgF, integrating it over R^ after adding the resulting equations 
together implies 



i.i(||A,E||i. + ||A,F||i. 
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(3.8) = / ■ ^ ^Aq{h{g)u) ■ A,E, 

where we used the vector analysis formula V • (f x g) = (V x f ) • g — (V x g) • f . 

Combining with the above identities (|3.7p - (|3.8p . with the aid of Cauchy-Schwartz 
inequality, we get 

i^(||A,e||i. + l|A,u||i. + i||A,E|li. + ^||A,F||i.) + i||A,u||i. 

< ^\\Vu\\L^i\\A,g\\l2 + \\Agu\\l.) + l^||V^?||L~||A,^?|U2||Agu|U2 
+ ||[u,A,]V£)||L2||Ag^)||i2 + ||[u,Ag] • Vu||i2||Agu||i2 

A,]Vg|U2||A,u|U2 + ^\\[g, Ag]divu||L2|lAqe|li2 

(3.9) +e||A,(u X F)|U2||A,u||i2 + i||A,(/^(f?)u)|U2||A,E|U2, t e [0,ro]. 

Next, we may neglect the effect of relaxation term i||Aqu||^2, since it is only re- 
sponsible for the large time behavior of solutions to p.2p - (|3T3l) . Dividing (|3.9p by 
{{WAqgWl, + II Aqull^, + i|| A,E||2, + i|| A,F||2,) + e}5 (e > a smaU quantity), we 
obtain 



2 dt I 



1 1 \ 1 ^ 

I A,g||i. + II A,u||i. + ^||A,E||i. + ^l|A,F||i.) + ej 



< ^l|Vu||i^(||A,e|U2 + ||A,u|U2) + I— I||v^?||i^||A,u|U2 
+ \\[u,Ag]\/g\\L2 + ||[u. A,] • Vu|U2 + l—^\\[g^Aq]\/g\\L2 

(3.10) +-^\\ I^'' A,]divu||i2 + £|| A,(u X F)||i2 + ^\\Agihig)u)\\L2 

for t S [0,To]. Integrating p.lOp with respect to the variable t, then taking e — >■ 0, 
and using the estimates of commutators and continuity for the composition in the 
stationary case, see [TU], we arrive at 

2?n|A«W-||Lr(L=) 

<C2«'^||A,M^olU^+C / c,{,mg,u,F)\\l. d, 
Jo 

(3.11) +C f 2'n|(Ve,Vu)||i^||(A,e,A,u)|U2d<^, 

Jo 

where ||cg(t)||^i < 1, for all i G [OjTo]. Next, summing up (13. lip on g > —1 gives 



(3.12) \mi^^^.^^^<C\\Wo\\Bi,+C \\W{;,)rs.J,, te[0,To]- 

■ J 

Then it follows from Remark 12.11 and Gronwall's inequality that 

(3.13) W&L^^iBl,). 

Furthermore, it is just a matter of using the equations p.2p and Proposition 12. 5[ we 
deduce that 

(3.14) WteL^^^iB^-,'). 
Hence, the proof of Proposition 13.11 is complete. □ 
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3.2. Uniform a priori estimate and global existence. In this section, our 
central task is to derive a crucial (uniform) a priori estimate, which enables us to 
achieve the global existence of classical solutions to 



Proposition 3.2. IfWe Ct{B^^i) n C^{BlJ^) is a solution of ^S^^-^S^^ for 
any T > and < r, e < 1. There exist some positive constants (5i,/ii and Ci 
independent of (r, e) such that for any T > 0, if 



(3.15) 

then 



\m 









) 


+ 











(3.16) <Ci||M^o||b|_,. 



Having Proposition 13.21 thanks to the standard continuation argument, we can 
extend the local- in-time solutions in Proposition 13. 1[ and achieve the global existence 
of classical solutions to the system (I3.2p - p.3p . here we omit details, see e.g. [TU]. 
It follows from Remark [^TTl Proposition 12.51 and the imbedding property i ^ 
that W e Ci([0,oo) X R^) solves dX^-dSS]). The choice of Si is sufficient to ensure 



^g+ip > 0. Then according to RemarklXTl we know (n, u, E, B) e C^{[0, oo) x R^^ 



is a solution of (|l.ip - (|1.3p with n > 0. Furthermore, we arrive at Theorem 11.21 

Actually, the proof of Proposition |321 is to capture the dissipation rates from con- 
tributions of {g, u, E, F) in turn by using the high- and low-frequency decomposition 
methods. To do this, we divide it into several lemmas. 



Lemma 3.3. If W e CriB^i) n C].{B^J^) is a solution of fSl)-fl3) for any 
T > and < t, e < 1, then the following estimate holds: 



(3.17) 



<C\\Wo\\b,,+C 



\w\\ 



L?,(i3|,) 



where /i2,C are some uniform positive constants independent of (T,e). 

Proof By integrating p. 91) with respect to t G [0,T], with the help of Cauchy- 
Schwartz inequality, we have 



L2 



1 



|A,E||i. 



1 



A,F||i. 



1, 



< 



^l|Vu||i2^(ioo)(||Age||i|(i2)||Aqg||ioo(i2) -I- ||AgU||i2^(i2)||AgU||i=o(i2 
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1, 



2 nV£'||L-(L~)l|A,£'||L2^(i2)||A,u||i2^(i2) + ||[U, Aq] • V^?||i2^(i2)||A5£»|li|(i2) 

7-1 

'II[U, Aq] • Vu||i,2,(i2)||AgU||i|(i2) + — A5]V6l||i2^(i2)||A,u||i2^(i,2) 

7-1 

— r— A,]divu||i2^(i2)||Ag£<||i2,(i2) +e||A,(u x F)||i2^(i2)|| A5u||i|(i2) 
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(3.18) 



|Ag(/l(£')u)||il^(i2)||AgE||ioo(i2). 



There exists a constant /i2 > independent of (r, e) after multiplying the factor 2'^'^"' 
on both sides of p.lSp . such that 



M2 



229'^||A„u|l? 



<C22'"^||A,Wo|li2 

+^{ll"IIZ|(sj,)(^'llfllZ|(i3j,)ll^'ll 



(3.19) +c^||uxF||~^^^^^^^|lu|l~^(^^^^^+c2||M^,)u||~,(^„^^^^^ 



where we used Remark 1 2. 11 Lemma [2.61 and the smallness of e(0 < £ < 1); Here and 
below C > denotes a uniform constant independent of (t, e); {cg} denotes some 
sequence which satisfies |l(cg)||;i < 1 although each {cq} is possibly different in p. 191) . 

Then, with aid of Young's inequality (-/Tff < (/ + 5)72, f,g> 0), it follows from 
Proposition 12.51 and Proposition 12.71 that 

M2, 



2^^\\AqWh^^L.) + ^^2«'^||A,u|L^(i.) 
<C2'"'\\AgWo\\L- 



1 



+CcqJ\\u\\r^,^^ , — ||u 
(3.20) +Ccy 



lEl 



^ll^?IIZ|(Bj,) + ;^ll"llz|(Bjj 



lL-(i3J,) 

Hence, summing up p.20p on (7 > — 1 yields 

/M2 



liw^ll 



< 



which is just the desired inequality (13. 17^ . □ 

Lemma 3.4. If W e Ct{B^.i) n Ct(^2,T^) a solution of figi-flj) /or anj/ 
r > anrf < T, £ < 1, then the following estimate holds: 



<C{\\{q,u)\\ 



+ ll(£'o,Uo)|lsj J 



(3.21) 



-C{;^||u||. 



ll(e,u,F) 



where C is a uniform positive constant independent of (r, e) 
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Proof. Set 



u y ' -J^ ^ \ '4'ej uUnxn 
{Inxn denotes the unit matrix of order N 

and ej is iV-dimensional vector where the jth component is one, others are zero). 

Then the first two equations of p.2p for Wj can be written as the following vector 
form 

(3.22) d,W, +J2A^iu)d.^Wj ^ _m:&,v7+G ) ' 

where G := -(E + eux (F + B)). 

To capture the dissipation rate of g, we make the best use of Shizuta-Kawashima 
skew-symmetric condition in Fourier spaces, which was developed for general hyper- 
bolic systems of balance laws [HI [29] . Thanks to the isentropic Euler equations p.22p , 
the concrete information of skew-symmetry matrix K(^) is well known {e.g. see [S]), 
which is very helpful to estimate the coupled electromagnetic field (E,B). Now we 
state the structural condition. 

Lemma 3.5 (Shizuta-Kawashima). For all ^ e R^, ^ ^ 0, there exists a real 
skew- symmetric smooth matrix K{^) which is defined in the unit sphere S^~^: 



(3.23) K{0 
such that 



0^ 



(3.24) K{i)Y^i,A\{^)-- ' 



-^f^ 



where Aj is the matrix appearing in the system \3.2't 
First, we rewrite p. 221) into the linearized form 



N , 2^ J- 

(3.25) dtWi^Y. 4 (0)^-. M^/ - e + f _ u _ v^^^^ G 

i=l \ r 2 

where 

(3.26) g = 5]{Aj(0)-A}(u)}9,^.l^,. 
Applying the operator A, to the system (I3.25P gives 

N 

(3.27) =A,g+f ^A^ul'^^^^'^^" 



2 



Ag(f?Ve) + A<,G 
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Then we perform the Fourier transform (in the space variable x) for (I3.27p . muhiply 
the resuhing equation by —iT{AgWi)* K {(,){* represents transpose and conjugator), 
and take the real part of each term in the equality. Using the expression (|3.23l) of the 
matrix if (^) we obtain 

, N 

rIm((A;W^,)*if(C) - aT^/) + r(A7fF,)*if (0 ( ^ ^.^(O)) A,Wi 

= rIm((A7l?^,)*if(0(A;^)) -Im((Z^)^A;u) +rIm((Z^)^A;G 

(3.28) +2_^rIm(Vi. ||(A,('^vu))) - I_IrIm(S^^(A,I^^.))) . 
The skew-symmetry of K{^) implies the relation 

(3.29) lm(^i\jWirKiOj^A^i'^ = ~lm(^i\jWir KiO\jWi 

Substituting p.24p into the second term on the left-hand side of p.28p . it is not 
difficult to get a lower bound. Indeed, we have 

, N 

(3.30) > ~lm(^iA:;WirKiOA^i) + ^miKiWi]^ - 2ViT|^||A;up. 

With the help of Young inequality and the uniform boundedness of the matrix K{^){^ ^ 
0), the right-side of p.28p can be estimated as 

Tlm((A7l?^,)*if(0(A;&)) -Im((Z^)|| A^u) +rIm((Z^)^A;G 

+ l_Irfni(Z^. ||(A,(^vu))) - I_IrIm(A^^(A,^^.))) 

< -f + T\[Q^,\\K;g\ -t- cr|A;u|i(A,(^vu))i 

(3.31) +CT|A7£-||(A,'(^Vf7))| +rIm((Z^)^A;G), 

where C > is a constant independent of (t, e). Combining the equality p.28p and 
the inequality p.30p - p.31l) . we deduce that 



<7(l^l + ^)|A,u|2 + r|A,iy,||A,c;| 

-f Cr| A^ull (A,(^vu)) l-f Cr| A^ell (A^T^Vp))! 
(3.32) +Tlm((A7g)^A;G) - ^-Iml^i^^Wiy K{0\;Wi 
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Multiplying (I3.32p by |^| and integrating it over [0,t] x R^, then using Plancherel's 
theorem yields 

^ \\A,VWi\\l.d, 
<- f {\\Agu\\l, + \\AgWu\\l,)dc^ + Ct f \\AgWWi\\L4\G\\L^d<; 

T Jo Jo 

+Ct ||A,Vu||i2 II A,(£-divu)||i2d^ + Cr / \\AgW g\\L2\\Ag{gW g)\\L2dc; 
Jo Jo 

(3.33) J (^iA^g)e^)d^d^^^lm J {^{(^(AgWi)* Ki^A^^Wi^d^ 

The matrix A'(^) is uniform bounded when f e R^C'C 7^ 0), thus we have 
-^Im J \^\(^{\WirK{0\iWi)d^\'^ 
<CT[J{l + \^f)\A,Wi{t)fd^ + J{l + \^f)\A,Wi{0)fd( 



(3.34) < Ci\\AgWiit)\\l2 + \\A,VWimh + \WWi{0)\\l. + || A,VW^/(0)||i2), 

where we used the smallness of r(0 < r < 1) in the last step. 

Next we turn to estimate the coupled electromagnetic field (E, B): 



r j^l^j (^{A^g)e \Gy^d<; 

j^^l^j ((A^g)^ ^)d^d^ 
-TS Im j ((A7e)e^A,(ux7F + B)))d^d^ 



(3.35) =Ii+l2, 

where the first term Ii is estimated as follows 



Im j [{A^g)C AqE>jdidq 
{- [ij {{\g)^^\^)d^d, + ^ J {{A;g)e^^d^d,} 



= -t{2t:)" j j AqVg ■ AqEdxdc, 
= t(2tt)^ J J AggAgdivEdxdi; 
= -t{2tt)^ J AggAgihig) - h{0))dxd^ 

= -T{Poj)--^n-^{27T)'' f ||A,e||i2d^-r(2^)^ / I AggAg{h{g)g)dxd<, 

Jo Jo J 
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(3.36) <-T(Po7)"^«^(27r)^ f \\Agg\\l,d^ + Cr T || A«(/»(e)e)IUHI A,^^|U2d^. 

Jo Jo 

Here, h{g) = h'{eg)d<; — {PQ"f)^in^^ is a smooth function on {g\-^^e-g + > 
0, e £ [0, 1]} satisfying h{0) = 0. 

In a similar way, I2 is estimated as 

I2 = -rej^ Im ^ ((A7e)C^A,(u xTf + B)))^^* 

(3.37) -'-^^^l (||A,div(uxB)|U2 + ||A,div(uxF)|U2)||A,^)|U2d<^. 
Thus, combining with p.33p - p.37p . we get 

^ f\\A,\/Wj\\l.d, + r{P,jr^n'^{2n)^ f \\A,g\\l,d<; 
^ Jo Jo 

< C{\\AgWi{t)\\l, + \\Aq\/Wi{t)\\l, + ||A,M^/(0)||2, + \\Aq\/Wi{0)\\l.) 
+ - [ (||A,u||i2 + ||A,Vu||i2)d^ + Cr / \\AgWWi\\L2\\Agg\\L2d<, 

T Jo Jo 

+Ct ||AgVu||L2||A,(edivu)||i2d^ + Cr / \\AqW gU-Al^qiQ^^ Q^lL^dc; 
Jo Jo 

+Ct f \\Ag{h{g)g)\\L2\\Agg\\L2d^ 
Jo 

(3.38) +CTe (||A,div(uxB)||i2 + ||A,div(uxF)||i2)||A,f?||i2d^. 



< C(22^||A,Ty,|li^(^2) + 2^''\\A,Wjml.) + c{ — \\A,ur 



Recalhng Lemma [2T2l we have 

||A,V/|U2 «2«||A,/|U2 (q>0). 
Note that this fact, from p.38p . we get the high-frequency part of || Agg||^2(^2)(q > 0): 

^2^^||A,,||i2(,2) 

■229, 

+re||A,e||i|(i2)||A,div(u x 6)11^2^(^2) + 2«r||A,l^,||i2^(i2)||Agg||^|(i2) 
+2V||A,u||i^(i2)||A,(gdivu)||^^(i2)+2«r||A5^.||i2^(^2)||A,(^)Vg)|L2^(i2) 

(3.39) +T\\A,Chig)g)\\L-^iL2)\Wg\\LUL2) + T£||A,div(u X r)||i2^(i2)||A,e||i^(i2)} 
and the corresponding low- frequency part: 

r(Po7)-^n'^(27r)^||A_ie||i2(^2) 
< C(||A_iiy,||i^(^2) + ||A_iW^,(0)||i2) 

-fCj^ll A_iu|||2^(^.) + A-if?||L|(L2)|| A_idiv(u X 6)11^2^(^2) 
-hT||A_iW>||i2^(i2)||A_ig||i|(i2) +r||A_iu||£2^(i2)||A_i(£»divu)||i|(i2) 

+ ^l|A-l6l||L2^(L2)||A_l((?V(?)||i2^(i2) +r||A_l(/l((3)e)||i|(i2)||A_l£»||i|(i2) 

(3.40) +Te||A_idiv(u x F)||i|(i2)||A_ie||^|(i2)}. 
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To conclude, we combine (I3.39|) - (l3.40p and multiply the factor 2^'^('^ on both sides 
of the resulting inequality to obtain 

t2o(T|| a _I|2 



<cci{\\wr 



2^ 



||div(uxB)||~^^^_,^ 



+TC': 



\\hig)g\\ 



(3.41) 



-e||div(u X F)|| 



lUb: 



where {cq} denotes some sequence which satisfies ||(cq)||/i < 1. 
By employing Young's inequality, we are led to the estimate 



(L2 



<CCgi\\Wl\\ 



L^(B- ) 



^(BJl)ll"llL|(BJi) 
1 

■ ull 



ll^llL-(i3J,)ll"llL|(SJJ 



(3.42) +c«y'll(^?,F)||j;^(^„^)(V^||^?||j;^(^.^) + -^||u|l~,^(^^^^^) }. 

In the end, with the help of the smallness of (t, e), summing up p.42p on g > — 1 
concludes the inequality (I3.2ip readily. □ 

Lemma 3.6. If W e CT(Sf^i) nC^(S^7^) is a solution of fS^j-fSJ) for any 
T > and < r, e < 1, then the following estimate holds: 



<C(||(u,E,F 



L^(Bl^) 



||(uo,Eo,Fo)||sjJ 





+ 


VF 




LliBl^) 





(3.43) +./||(^>,u,F)|| 



VF 



Ll(Bl-' 



where C > is a uniform constant independent of (r, e) . 

Proof. A nice "div-curl" construction of Maxwell's equations of p.2p enables 
us to obtain the high-frequency part of E. Indeed, by applying Aq to both side of 
V • E = —h{g){q > 0), integrating it over after multiplying V • A^E, in virtue of 
Holder's inequality, we obtain 

l|V-A,E||i. 

(3.44) < C{{Po^)-in-^\\Aqg\\L2 + \\Aq{h{g)g)\\L2}\\V ■ A^EjU^, 
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where the function h{g) is defined by (j3.36p . 

On the other hand, applying Aq{q > 0) to the fourth equation of p.2p and 
muhiplying the resuUing equation by V x AgE, after integration by parts, yields 

||V X A^Elli^ J dtA.F ■ (V X A,E) 

(3.45) = e j{V X atA.F) • A,E. 

Substituting the third equation of p.2p into p.45p . by Cauchy-Schwartz inequality, 
leads to 

||V X A^Elll^ + ||V X AgF||2, 
<e— /(V X A,F) • AgE + n£||A,u|U2||A,(V x F)\\l2 

(3.46) +e\\A,{h{g)u)\\L2\\A,(y x F)\\l2. 

Combining with p.44p and (|3.46p . it follows from the elementary relation 

||Vf|U2«||V-f|U2+||Vxf|U2 



that 



|VA,E|12, 

(Pi 
d 



< C{{Po^)---n-^\\Agg\\L2 + \\Agihig)g)\\L2}\\V ■ A,E||i2 



dt 



■ y"(V X A,F) • A,E + ne\\Agu\\L2\\Aq{^ x F)\\l2 

(3.47) +e\\Agihig)u)U2\\A,{V x F)|U2. 
Note that g > 0, from Lemma [521 we further get 

re2^^\\A,F\\l2 

< CT\\A,g\\l2 + CTs\\A,{h{g)g)\\L22'^\\A,E\\L2 

+Te2A y"(v X A,F) • A,E + Ct£2||A,u|U2||A,(V x F)|U2 

(3.48) +T£2||A,(/i(e)u)|U2||A,(V X F)||i2, 

where the smallness of e is used. Integrating p.48p in i G [0,T] implies 

r£22^|iA,E||i2(i2) 

< re229||A,F||io,(i2)|lA,E|Uo,(i2) +r£22«||A,Fo||L2|lA,Eo|lL2 
+Cr|| A,e||i^(^.) + CTe||A,(Me)g)|L^(i2)2'|| A,E|L^(i2) 
+Cr£2||A,u||i^(i2)||A,(V X F)||i|(i2) 

(3.49) +Te2||A,(/i(^.)u)|L^(i2)||A,(V x F)|L^(i2). 

On the other hand, the desired low-frequency of E can be deduced from the 
Lorentz field in the Euler equations of p.2p . Using the second equation of p.2p . we 
have 

(3.50) ut + E = -tpVg - - - u • Vu - ^^^fiiVp - eu x (F + B). 

T 2 
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Applying the operator A_i to (|3.50l) iraplies 

atA_iu + A_iE = -i/^A-iVp - ^-^^ - A_i(u • Vu) 

T 

(3.51) -^^A_i(f?Ve)-£A_i(ux (F + B)). 

Multiplying (I3.5ip by TeA_iE and integrating the resulting equation over R^, we 
get 

TS— j A_iu- A_iE + T£||A_iE|||2 +Te7/;(^'o7)"^^^'^l|A-if?||i2 
— nre ||A_iu|||2 +r j A_iu- A_i(V X F) + re j A_iu • l\q{h{g)u) 
+Te'4) j A_i(/i(£»)^?)A_i£»-e j A_iu • A_iE 
-T£^ j (A_iu X B) • A_iE - T^y" A_i(u • Vu) • A_iE 

(3.52) -^-y-'^'^ j A-i(£'Ve) • A_iE - re^ j A_i(u x F) • A_iE, 

where we have used the third equation of p.2p . From Cauchy-Schwartz and Young's 
inequalities, we arrive at 

j A-iu-A_iE+^|lA_iE|li. 

< c(i||A_iu||i. +t||A_iu|U2||A_i(V x F)!^^) 

+Cre(||A_l(/^(^))u)|U2||A_lu|U2 + ||A_l(/l(^?)^))||i2|| A_if?|U2 
+ ||A_i(u- Vu)||L2||A_iE||i2 + ||A_l(gV^^)||i2||A_lE|U2 

(3.53) +e||A_i(ux F)|U2||A_iE|U2' 

Then integrating p.53p in t e [0, T] gives 

^£||A-lE||i2(^2) 

< Cre(||A_iu|Uoo(i2)|lA_iE|Uoc(i2) + |1 A_iUo||l2||A_iEo|U2) 
+c(^||A_iu||2^(^,) +r||A_iu||i2^(i2)||A_i(V x F)||i2^(i2 

+CTe(\\ A_l(ft.(£))u)||i2^(i2)||A_lU||i2^(i2) + ||A_l(/l(£l)£))||i|(i2) 

■\\^-iQ\\lI(l^) + l|A-i(u • Vu)||i|(i2)||A_iE||i2,(i2) 

(3.54) +l|A_l(^.V^.)|L^(i2)||A_lE|L^(i2) +e||A_i(u x F)||i^(i2)|| A_iE||i^(i2) 

Therefore, by combining with the high-frequency estimate p.49p and low-frequency 
estimate p.54p . we infer that for g > — 1, 

Te22«||A,E||i2(^2) 
< C229(||(A,u, A,E, A,F)||2^(^2) + ||(A,uo, A.Eq, A,Fo)||i2) 
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+Cts{\\ 

+ (||A,(u • Vu)||i2^(i2) + ||A,(^>V^?)||i2^(i2) + 2''\\Ag(h{g)g)\\L2^(^L-^ 
+e||A,(uxF)||i^(i2))||A,E||i|(i2) 

(3.55) +T£2||A,(/.(f,)u)L|(i2)||A,(V X F)\\LliL-)}- 

By multiplying the factor 2^'?'°'^^^ on both sides of p.55p . we gather 
r£22«'^||A,E||i.(^.) 

X Fllr. 

(3.56) +£||u X F\\~,^^^.^^^.^)mi2^^^^^^ X F||~.(^„_.,}, 

where {cq} denotes some sequence which satisfies ||(cq)||;i < 1. 
Then it follows from Young's inequality that 

V^2*n|A,E|U^(i.) 
< Cc,(||(u,E,F)||~^^^„^^ + ||(uo,Eo,Fo)||bjJ 

(3.57) +c,^^^^^(-L||u||~, + 7^11^ ^ ^llz^(Bn^))' 

where we has used Proposition 12.51 Finally, we sum up p. 571) on q > — 1 and deduce 
the inequality p.43p immediately. □ 

Lemma 3.7. If W e Ct{B^.i) n Ct(^2,T^) « solution of ([13)-([2g) /or any 
T > and < t, £ < 1 , then the following estimate holds: 

7|II^F|lzi(B,7^) 
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<C(||(E,F)|| 



(3.58) +C{ 



L~(B| ) 



||(Eo,Fo)||b.J 



Ikll 



VF 



where C > is a uniform positive constant independent of (t, e). 

Proof Multiply the third equation of p.2p by — Aq(V x F)((7 > —1) and integrate 
the resulting quality over R^. Then integration by parts implies 



2 

L2 



(3.59) 



- / A,(V X E) • A,F + -||A,(V X F 
j A,(V X E) • A,Ft - n j A,u • A,(V x F 

- j A,{h{g)u) ■ Ag(V X F). 



into the first term of p.59p . by Cauchy- 



Substituting the fourth equation of 
Schwartz inequality, leads to 

d_ 

di , 

(3.60) < n||A,u|U2||A,(V X F)|U2 + || A,(M£))u)|U2 1| A,(V x F)||i2. 

Due to the incompressible condition of F in (|3.2p , by integrating p.60p with respect 
to t S [0,T], we easily derive 



A,(V X E) • A,F + i||A,(V X F)||2, + i||A,(V x 



1, 



A„VF| 



< l|A,(V X E)|U=e(i2)||A,F|U^(i2) + ||A,(V X Eo)|U2 1| A,Fo|U2 

+ "l|AgU||i2,(i2)||A,(V X F)||i2^(i2) 

(3.61) +|lA,(/^(^))u)L^(i2)||A,(V x F)L2^(i2). 

Noticing that the regularity of E in the assumption of Lemma l377l we multiply p.6ip 
by the factor 229('^-i) to get 

l229('T-l)||A VF||2 



(3.62) 



^4\\hig)u\ 



where {c^} denotes some sequence which satisfies |i(cq)||/i < 1. 

Furthermore, we apply Young's inequality to p.62p and obtain 



i=2^(--i)||A,VFL^(^2) 
<Cc,{||(E,F)||~^^^.^^ + ||(Eo,Fo)|1bjJ 
'{c,^\\ 



-c 



(3.63) 



^Cq^\\h{g)u\ 



IV X 
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Finally, after summing up (|3.63p on g > — 1, it follows from Proposition 12.51 that 
1 



<C(||(E,F)|| 



|(Eo,Fo)||sjJ + c{^||u|| 



This is just the inequality p.58p . Hence, the proof of Lemma [H771 is complete. □ 

Remark 3.2. In the proof of Lemma [5771 the dissipation rate of F is not available 
due to the absence of low-frequency estimate on || A_iF||2:,2 (^2). This is a key reason 
that Chemin-Lerner's spaces with critical regularity are first introduced to establish 
the global existence of uniform classical solutions. Otherwise, we need to add a 
little regularity in order to ensure that the Besov spaces (in x) are still continuously 
embedded in C\li'^) spaces. For the similar details, the reader is referred to j9]. 
Having these lemmas proved above, the proof of Proposition 13.21 can be finished. 



Proof of Proposition [H Combing ([3T7| . ([33T|) . (|03l) and ([3381) . we end up with 



l|w^llz^(^.^)+i^l^/^ll^^llz.(BJJ + v/vll"llI|(Bf,,) 



+is:2V^||E|| 



\{g, u,F) 



+CA'2{(||(u,E,F)||~^^^^,^^ + ||(uo,Eo,Fo)||bjJ 

1 

(e,u,F)||~^^P„ J (V^e, ^u, 



^VF 



(3.64) +^I|VF||~2 

where K\^Ki and i^a are some uniform positive constants (independent of (t, e)) to 
be determined. In order to eliminate the terms ||(£i, u, E, F)|j~^^^^ j, ||\/r£'ll2^2 (-^.r 

\\vil ypT\\j.2 , , and ||VF/ye||7-2 /R,^-iN arising in the right-hand side of p.64p . we 

2.\) t'- 2,1 1 
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may confine the constants to the fohowing cases: 
Furthermore, it is not difficult to obtain 



+if2\/^||E|| 



2 4^11 "i?^(Sf,i) 

^ IVFII 



+CK^{\\i go, uo)||bj^ + ^Wig, u, F)\\~^^^^^^ 



(^?,u,F)|| 



y/rg, —;=u, y/rsE 

\ T 



VF 









}+CX3{||(Eo,Fo)||bj_^ 



;,)(;^ll"llz^(sj,,) + 7fll^F|IZ|(Bn^))} 



< 



VF -1 



(3.65) < CIIW^oIIbj, +Cv^{||(V^^.,^u,V7iE 



) 


+ 


VF 









where we used the a priori assumption p.lSp in the last step of p.65p . 
Finally, we choose the positive constant 5i such that 



1 < mm < , , iv 2 , > , 

I 2 4 2 i 



then the inequality p.l6p follows immediately. This finishes the proof of Proposition 
eventually. 



3.3. Non-relativistic limit. In this section, we justify the non-relativistic limit 
of the system (fLTjl - lfO)) with t 1. 



Proof of Theorem\1.4\ For any fixed T > 0, let (n^, u"^, E'^, B"^) be the global solution 
of (|l.ip - (|1.3p given by Theorem ll.2l It follows from the uniform energy estimate (|l.lip 
and Remark |2. II that 



(3.66) 
(3.67) 
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(3.68) B'-B eL^{Bl,), ^ lUb^j'), 

uniformly in e. Note that p.68p . we deduce 



(3.69) <CVe-^0, as e ^ 0. 
That is, 

(3.70) {VB^}->0 strongly in L^B^^), as e ^ 0. 
Moreover, with the help of we have 

(3.71) (nlul) e L^B^'). 

(3.72) ViE? e 4(S2,l'), 



uniformly in e. 

According to p.66l) - p.68l) and p.7ip - p.72p . it can be derived from Proposi- 
tion 12.31 and Aubin-Lions compactness lemma in [261 that there exists some func- 
tion (n",u",E°) e C([0,oo),n + B£j xC([0,oo),B£i) xC([0,oo),B^_i) such that the 
sequences (up to subsequences) as e — > 0, it holds that 



(3.73) Kl^n" strongly in C{[0,T], {B^%oc), 

(3.74) {u^}^u° strongly in Ci[0,Tl {B^%oc), 

(3.75) {V^E'j^O strongly in C{[0,T], {B^f)i,,), 

(3.76) {E^}^E° weakly* in L^iB^-^), 

(3.77) {B^I^B weakly* in i??(BJi), 



for any T > and S G (0, 1). Thus, in the system (fLT|) - ([L3l) . the unifo rm bounded 
properties p.66p - p.68p as well as the convergence properties p.69p and p.73p - p.77p 
allow us to pass to the limit e — in the sense of distributions, which implies that 
E") is a global weak solution to the Euler-Poisson equations (11.61) satisfying 
(|1.12p . This completes the proof of Theorem 11.41 

3.4. Relaxation limit. In this section, we prove the relaxation limit of (|l.ip - 
(fOl) with e = 1. 

Proof of Theorem \1.5i From the scaled variable transform ()1.7p and the uniform en- 
ergy estimate (jl.lip in Theorem ll.2[ it is shown that (n'^, u^, E^, B^) is a unique 
global solution of the system (|1.8p and (|1.13p . furthermore, for any fixed T > 0, we 
have 

(3.78) n^-neL^(S^_i)nL^(B^_i), 
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(3.79) ru^ e L^{Bl,), e L^ABli), 

(3.80) eL^iBl,)nLUBl,), 

(3.81) B--BeL??(B2%), ^ei|(S2"l'), 

uniformly in r. Similar to p.69p . it follows from p. 811) that 

(3.82) {VB^}->0 strongly in L^{B^^), as r ^ 0. 
Moreover, from the equations (|1.8p . we conclude that 

(3.83) < € iK^a"!'), 

(3.84) r2u[ G ii.(i?27i), 

(3.85) v^E[ e £|(S2"T')> 
uniformly in r. 



Together with ((3J8l) - ((3:8T]) and jSlSll-dSH]), it follows from Proposition [23] and 
Aubin-Lions compactness lemma in 26 that there exists some function {Af,U,£) G 
C([0, oo),n + Sf i) X L2([0, oo),Bf i) x C([0, oo), j) such that the sequences (up to 
subsequences) as r — >■ 0, it holds that 



(3.86) {n^j^Af strongly in Ci[0,T], {B^%,,), 

(3.87) {r^u'l^O strongly in Ci[0,T], {B^-,%oc), 

(3.88) {u"'}^^^ weakly in L^Bl^), 

(3.89) {V^E^j^O strongly in C{[0,T], {B:^-^%,), 

(3.90) {E^j-^S weakly* in L^{Bl{), 

(3.91) {B^I^B weakly* in Lt(^2.i), 



for any T > and 5 G (0, 1). Thus, the uniform bounded properties (|3.78p - p.8ip 
as well as the convergence properties p.82p and p.86p - p.9ip allow us to pass to the 
limit T — 7> in the system (jl.Sp and ()1.13p in the sense of distributions, which implies 
that (A/", £) is a global weak solution to the drift-diffusion equations (|1.9p satisfying 
(|1.14p . Hence, the proof of Theorem 11.51 is complete. 
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3.5. Combined non-relativistic and relaxation limits. In the last section, 
we perform the combined relativistic and relaxation limits of (jl.ip - (ll.3p . 

Proof of Theorem \1.6l Combined with the scaled variable transform (I1.7P where the 
superscript r is replaced by (r, e) and the uniform energy estimate (jl.lip in Theorem 
11.21 it is obtained that (n^^''^^ u^'^-^), E*-'^''^-', B^'^''^'') is a unique solution of the system 
(jl.lOP and (I1.13P with the superscript r replaced by (r, e) , furthermore, for any fixed 



T > 0, we infer that 

(3.92) -ne i?? (^2,1) <^ ^t(52,i), 

(3.93) ru(^^^) e (^2,1), u^^-') e ^1(^2,1), 

(3.94) E^^'^) e L^{Bli), ViE^^'") e L^ABli), 

(3.95) B(-'^)-Bei??(i?2^,i), e lUb!^-'), 

uniformly in (t, e). The relation p. 951) turns out to yield 

(3.96) {VB("''^)}^0 strongly in L^B^^), asT,e^O. 
Moreover, using the equations (jl.lOl) . we get 

(3.97) n^-'^) e L^B^,-'), 

(3.98) T^n^r'' e L^BI^ri'), 

(3.99) V^E^^) e L^B^J'), 
uniformly in (t, e). 



As previously, it follows from the standard weak convergence methods and com- 
pactness lemma in [55] that there exists some function {J\f,U,£) e C([0, oo),n + 
^2,1) X -^^([Oj 00)7-^2.1) ^ C([Oi 00): ^2,1) such that the sequences (up to subsequences) 
as T, £ — > 0, it holds that 



(3.100) 




strongly in 


c([o,r],(i?2-7*)ioc), 


(3.101) 




strongly in 


C{[0,T],{B^')i,,), 


(3.102) 


{uKs)}^Z^ weakly 


in LUBI^), 


(3.103) 


{^/7iE("-')} ^ 


strongly in 


C{[0,T],{B^%,,), 


(3.104) 




£ weakly* 


in L^{Bl,), 


(3.105) 


{B(--^)} - 


B weakly"^ 


in L??(i32%), 



for any T > and (5 G (0,1)- Thus, in the system (jl.lOl) and (jl.lSp . the uniform 
bounded properties p.92p - p.95p as well as the convergence properties p.96p and 
p.lOOp - p.lOSp allow us to pass to the limits r, e — >■ in the sense of distributions, 
which implies that {J\f, £) is a global weak solution to the drift-diffusion equations 
(|1.9p satisfying (|1.15p . This concludes the proof of Theorem 11.61 
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